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We summarize results presented at this conference with special empha- 
' sis on hard processes with jets and heavy quarks, soft particle production, 

\Q . small x structure functions and diffraction as well as heavy ion collisions 

and quark gluon plasma. 



■^j - ■ 1. Introduction 

o 

In high energy collisions of leptons, hadrons and nuclei we observe the 
production of many particles, mainly hadrons. The ultimate goal in mul- 
tiparticle production studies is the explanation of the hadronic phenomena 
within QCD, the basic theory of the strong interactions, along with the 
other interactions of the standard model and possibly beyond. 

A basic problem in the QCD study of multiparticle production is the 
matching of parton and hadron dynamics relevant in the respective weak 
and strong coupling regimes of QCD. A class of inclusive observables in 
hard collisions can be computed perturbatively in terms of the running 
coupling constant a s thanks to the celebrated asymptotic freedom pQ. In 
practice, it is often required to include some additional non-perturbative 
input from other sources {e.g. Parton Distribution Functions). Hadronic 
phenomena can be systematically analysed within lattice gauge theory in 
sufficiently simple problems. The description of genuine multi-hadron pro- 
duction requires, in addition to the hard QCD part (shower calculus), some 
phenomenological approaches: specific hadronization models or such simple 
ideas as "parton hadron duality". At present, the systematic approach to 
multi-hadron production based entirely on QCD remains a dream. Is it just 
a problem of complexity or do we need a fundamentally new insight? 

In this talk we concentrate on the four topics mentioned in the abstract 
which represent different variants of the interplay of hard and soft interac- 
tions and we emphazise shortly other topics. I am sorry for the incomplete 
coverage of the many interesting presentations at this conference. 
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2. Hard processes with jets and heavy quarks 

Results are reported from TEVATRON, HERA, RHIC and LEP accel- 
erators. Of central importance is the comparison with fixed order pertur- 
bation theory to test the universality of the coupling constant a s (Q 2 ) in 
all processes. Here the goal is to improve the accuracy of the calculation, 
in particular by the extension beyond Next-to-Leading-Order accuracy, and 
to compute within the QCD framework observables of higher complexity. 
An important goal is also the discovery of new physics either through the 
deviation of experimental results from the precision calculations or through 
better understanding of background processes. 

Top quark production (D. Bauer) has been studied at the TEVATRON 
pp collider where in the new RUN II the cms energy is increased from 1.8 
to 1.96 TeV. Top quarks at these energies are produced primarily in pairs 
and they decay through t — > Wb. The various decay channels for the W 
have been analysed by both experiments and consistent results have been 
obtained in eight channels, all compatible with the theoretical computations 
beyond NLO of a top = 6.77 ± 0.42 pb at m t = 175 GeV |2j. A new value 
for the top quark mass has been presented by DO [3| from a reanalysis of 
their earlier RUN I data: m t = 180.1 ± 3.6(stat.) ± 3.9(syst.) GeV with 
considerably reduced errors and this result increases the world average by 4 
GeV to mt = 178 ± 4.3 GeV. The best RUN II CDF result so far obtained 
is m t = 177.8^5'q ± 6.2 GeV. No single top production has been observed 
and a t < 8.5 pb at 95% CL (RUN II, CDF). 

Single top production (S.D. Ellis). The observation of this process would 
determine the CKM matrix element Vtb and be important in other searches 
(Higgs, new particles, including extra scalar bosons or gauge bosons, new 
quarks). The expected rates are still below the presently achieved sensi- 
tivity. Therefore some strategies to obtain an improved signal/background 
ratio are proposed, especially by using the "signed rapidity" variable which 
takes into account the fact that processes with VF's are not separately C or 
P invariant. 

Jet Production at the TEVATRON (L. Sawyer) is measured now at the 
higher cms energy by CDF and DO. Single inclusive jets are studied up to 
transverse momenta pr ~ 550 GeV, that is 150 GeV higher than in RUN I 
and di-jet masses up to Mij ~ 1400 GeV, also the azimuthal decorrelation 
of the two jet events which is an effect of 0(a s ) 3 is measured. These results 
are found in an overall consistency with the NLO QCD predictions [IJ |5] 
in the extended regime of energy and transverse momentum within the sys- 
tematic errors which are dominated by the energy scale uncertainty. Only in 
extreme kinematic regions of the azimuthal decorrelation deviations appear. 
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Fig. 1. Distribution of the jet for the jet of lowest energy in the W+ > n jet 
sample, for n=l (top) to n=4 (down), in comparison with LO QCD calculations. 

Production of W,Z,j + jets (A. Cruz) allows further tests of pQCD 
at large momentum transfers using different techniques. The inclusive W 
production cross section rises from 2.38 ± 0.24 nb at 1.8 TeV to 2.64 ± 0.18 
nb at 1.96 TeV and is well described by the NNLO QCD results of 2.5 
and 2.73 nb respectively. For the production of W together with n jets the 
distribution of the jet transverse energies have been measured. There is 
good agreement with a calculation (ALP GEN jSJ) based on leading order 
QCD Matrix Element combined with parton shower development (HERWIG 
[7j), see Fig. ^ Here the experimental uncertainties from the jet energy 
scale are comparable in size to the theoretical scale uncertainties. Central 
di-photon production has been measured up to masses of about 30 GeV. 
The NLO QCD calculations (DIPHOX 8 ) involving the qq annihilation 
and quark loop diagram for gg — > 77 describe the data well in absolute 
normalization. The associated production of photon with heavy flavours 
(c,b) has been obtained as well for photon energies up to 60 GeV. It agrees 
with a LO calculation (Pythia 9 ). 

In summary for TEVATRON, no serious disagreement of the new data 
with QCD expectations at different levels of accuracy for the observables 
of different levels of complexity are met and so, at the same time, we do 
not obtain any signal of a new physics. The other problems discussed here 
concern the role of the "underlying event" {R.D. Field, see below) and the 
influence of the selection of jets by a particular algorithm ("cone" vs."-fT;r") 
(Andrieu). 

Monte Carlo generators, (L. Lonnblad) are being developed to generate 
parton final states beyond LO for present accelerators but also for the LHC. 
For problems like production of multi-jet or W+jet events a matching of 
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the NLO matrix element and parton shower has been achieved; also double 
scattering effects are being considered. 

Fragmentation functions (S. Kretzer) Dq^(z,Q 2 ) are compared in pp 
and e + e~ collisions as a tests the collinear factorization approach of QCD 
at NLO together with universality. New data obtained at RHIC on pp — * ttX 
are well predicted using the previous results ^3^D from e + e~ collisions on 
at the factorization scale and subsequent DGLAP evolution. 

Jet production at HERA (C. Glasman) has now been measured down 
to x ~ 10 -4 for momentum transfers Q 2 larger than a few GeV 2 . The 
Q 2 dependence of of two- and three-jet cross sections in NC DIS up to 
Q 2 ~ 5000 GeV 2 , obtained by ZEUS, is found in good agreement with 
NLO (0(a 2 ) and G(o? s )) predictions. The ratio of both results provides 
an accurate test of the theory and in particular a precise determination of 
the coupling constant a s (M z ) = 0.1179 ± 0.0013 (stat.)iaoM6 ( ex P-) which 
by itself compares well with the current world average a s (Mz) = 0.1182 ± 
0.0027 [12] but there is still a large theoretical error of (+0.0061, -0.0047). 

Of special interest, also in view of further applications to small x physics, 
is a test of the validty of the DGLAP approximation to the Q 2 evolution of 
structure functions. If jets are produced "forward" (in direction of incoming 
proton) and have transverse energies (E^ ) 2 of order Q 2 then the kinematic 
configuration is not in favour of intermediate gluon emission with strong hp 
ordering as is typical for DGLAP evolution and deviations from DGLAP 
based predictions are expected [T3|- Indeed, if jets are selected by the HI 
collaboration with 0.5 < (E J T et ) 2 /Q 2 < 5 and x < 0.004, then the "direct" 
NLO QCD predictions (DISENT |T1|), i.e. without photon structure, are too 
low by about a factor 3 for very small x < 0.001. Considerable improvement, 
although not full agreement, is obtained if a photon structure is included in 
the DGLAP calculation. If one of the two scales Et or Q 2 is large compared 
to the other one the DGLAP factorization approach works best. Whereas 
work continues to develop the approximations beyond DGLAP, there is no 
serious conflict with pQCD at the fundamental level. 

Spin physics. The spin program at RHIC (E. Sichtermann) aims at a 
measurement of of hard and soft processes with polarised protons. The 
single transverse spin asymmetry An of forward 7r° production has been 
observed at FNAL at yfs = 20 and the first results from RHIC show that it 
persists at yfs = 200 GeV. New results are also presented from the HERMES 
experiment at DESY (/. Gregor) on the study of the transverse single-spin 
asymmetries in semi-inclusive pion production in DIS. This allows the de- 
termination of the protons "transversity" distribution, which represents the 
degree to which the quarks are polarised along the proton spin transversely 
polarized to the virtual photon. 
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3. Soft particle production in hard collisions 

The classical applications of pQCD concern observables for hard pro- 
cesses where the hadrons are either summed over or collected into jets so 
that fixed order perturbation theory can be applied, as in the previous sec- 
tion. A further development, surprisingly successful, is the application of 
pQCD to observables calculated from the individual momenta of hadrons 
in the final state directly and in general involves a resummation of the 
perturbation theory. One has to include some assumptions on the tran- 
sition from partons to hadrons and eventually on some non-perturbative 
aspects of colour confinement. There are final state parton observables, like 
event shapes, which are infrared and collinear safe, i.e. their values do not 
change if a collinear or soft gluon is added. They are less sensitive to soft 
hadronization effects. More sensitive are multiplicity observables, such as 
particle flows inside jets or between jets which are not infrared safe. These 
observables are a testing ground for "parton hadron duality" ideas (review 
by Yu. Dokshitzer). 



3.1. Infrared and collinear safe observables 

Event shapes (Yu. Dokshitzer). These observables describe global prop- 
erties of hadronic final states, in e + e~ annihilation, for example, one defines 
"thrust", "jet mass", "broadening" and others. The analysis of these observ- 
ables in perturbation theory (an asymptotic expansion) leads to a descrip- 
tion which combines perturbative and non-perturbative aspects, the latter 
ones represented by a power correction oc (1/Q) |15| . This term involves 
an integral over gluon emissions at small transverse momenta k? where the 
coupling constant a s (kj,) is ill defined. It is assumed that this integral is fi- 
nite and - as the coupling itself - universal for the different observables. One 
therefore introduces the parameter ao = Jq 1 dkrOi s (kj^) at the matching 
scale m to describe the influence from the soft region [TJj|. This result, 
strictly obtained for partons, is then applied to the experimental hadronic 
observables assuming a duality between both descriptions. The calculation 
has been extended to the differential distribution of shape observables where 
the non-perturbative effects can shift or squeeze the perturbative spectra by 
an amount given by 1/Q. By now, the fits to event shapes in e + e~ with 
two parameters provide a competitive determination of a s with an error of 
~ 8% and the universality of the non-perturbative parameter ao ~ 0.5 at 
[ii = 2 GeV is confirmed within ~ 15%. 

Angularities (G.Sterman): A new class of event shapes allows specific 
tests for the non-perturbative corrections. From the angles Q{ of the particles 
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to the thrust axis and energies Ei one constructs the quantity ^2] 

Ta = T> E ^( sin ^)"(l " I cos ^l) 1 -" (1) 

w t 

which interpolates between thrust (a = 0) and broadening (a = 1). Again, 
one can separate a contribution from soft gluon emission which is then rep- 
resented by a non-perturbative "shape function" f a ,NP- It represents cor- 
rections of all higher orders in A/Q which should be more appropriate near 
the collinear limit. This represents a generalization of the correction ao/Q 
above, which corresponds to a shift of the distribution. Once determined at 
a particular energy one obtains predictions for other energies. 

Interjet radiation: non-global log's. In DIS and hadron-hadron collisions 
one is led to consider gluon radiation in part of the phase space exclud- 
ing a region around the beam direction, so energy flow or event shapes are 
non-global. Such observables obtain contributions from multi-soft emissions 
which lead to "non-global log's", single logarithmically enhanced contribu- 
tions [18j . The difficulty is that the number of jets is not fixed. This problem 
can be tamed by construction of a correlation of the energy flow with an 
event shape which fixes the number of jets |17l I19j . 

Automated resummation (G. Zanderighi): In order to facilitate the cal- 
culation of new shape observables, in particular for pp and DIS, a pro- 
gram (CAESAR |20j ) has been developed which, for a class of observables, 
yields resummed results in NLL order by combining analytical and numer- 
ical methods. The limitations concern certain properties in the infrared 
and collinear limit of the emission ("recursively IRC-safe") and the depen- 
dence on transverse momentum ("continuously global"). Old results have 
been reproduced, new observables derived, for example, the distribution of 
global transverse thrust Tj_ = ^ max^ T £V \p ti fiT\ constructed from the 
transverse momenta with respect to the beam axis. Such calculations open 
up the possibility to considerably extend the kinematic range of these QCD 
studies towards the higher energies of TEVATRON and LHC. 

3.2. Multiplicities, particle flows 

These observables are not infrared safe: emission of a soft gluon would 
increase the multiplicity by +1. Finite perturbative results for the parton 
cascade can be obtained by introducing a cut-off hp > Qo- For Qo 3> A the 
partons represent jets and Qo can be viewed as jet resolution in the sense 
of the "/c^-algoritm" , one can also take the small Qo ^ A and compare 
the resulting cascade directly with the hadronic final state in the sense of 
a duality picture ("Local Parton Hadron Duality" [21\). then Qo is a non- 
perturbative parameter. 
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Multiplicities in quark and gluonjets at LEP (K. Hamacher) and TEVA- 
TRON (A. Pronko): At LEP the multiplicity of gluon jets is determined 
from 3 jet events after subtraction of 2 jet events at a reduced scale (DEL- 
PHI) or from 3 jet events using a boost algorithm (OPAL [22] )• The <l uar k 
jet multiplicity is found directly from the total e + e~ multiplicity. Theo- 
retical results are obtained from resummed perturbation theory. One ap- 
proach is based on coupled evolution equations of quark and gluon jets in 
Modified Leading Logarithmic Approximation |21| which includes fully the 
y^-correction up to NLL order. These calculations reproduce the mul- 
tiplicity rise with energy. The ratio r = N g /N q obtains large corrections 
beyond MLLA and is reduced from the asymptotic value r = Ca/Cf = 9/4 
to r = 1.7 in 3NLLO |23| and to r = 1.5, observed at LEP, in the numeric 
solution [21] where the only parameters A, Qq are fit by the total e + e~ mul- 
tiplicity. Another calculation is based on the colour dipole model which 
treats the evolution of dipoles in NLL approximation and includes recoil 
effects [2H|- It describes the data well using an additional non-perturbative 
parameter. 

The CDF Collaboration has separated quark and gluon jets by analysing 
di-jet and 7 + jet events with known jet compositions. The multiplicity data 
are found generally well consistent with e + e~ results. The multiplicities 
reach the higher energies Q ~ 300 GeV where they also follow the 3NLLA 
expectations. 

Particles with low momenta in jets (A. Pronko), measured by CDF, 
show the so-called hump-backed plateau in the variable £ = log(l/x) [2J 
with the suppression of the low energy, large £ particles because of soft 
gluon coherence. Results on the ratio r(£) of these spectra for gluon over 
quark jets for £ > 3 approach ratios of r(£) ~ 1.8 =t 0.2 again in good 
agreement with OPAL results. Although larger than for the full jet result 
(r ~ 1.5) it is still below r(£) = Ca/Cf expected in |26| for this limit from 
the dominance of the primary gluon emission. This discrepency is likely 
due to the difficulty to obtain "pure" gluon jets, rather, soft particles are 
emitted from all participating jets of the event. This problem is avoided in 

Soft particle emission in 3- jet events in e + e~ (Hamacher). The particle 
multiplicity N$ in a cone perpendicular to the production plane is studied 
as function of the inter-jet angles Qij. The gluon radiation into this cone 
coherently emitted from the qqg "antenna" , normalized by a corresponding 
multiplicity N2 in 2-jet events, is given by the simple expression [2H] 



Ns = Ca = ICa 
N 2 C F n 4C F 



(l-COSOgg) + (1-COSOgg) ~ — 3 (l~COS Qgg) 

1\ C 



(2) 



The first two leading terms represent the dipoles along the qg directions, 
in close analogy to QED electric dipoles, except for the colour factors. The 



Fig. 2. Multiplicity ratio N 3 /N 2 in cones of 30° opening angle as function of r t in 
Eq.(j2J: a) for different inter-jet angles #3; b) averaged over #3, the dashed line is 
the expectation Eq.© with slope C A /C F , the full line is a fit (DELPHI |27|L 

formula interpolates for aligned partons between a colour triplet antenna 
(q against qg) and a colour octet antenna (qq against g) with the intensity 
higher by C A /C F . 

DELPHI has measured the ratio N3/N2 against the variable r t (see 
Fig. |2J) which shows the scaling behaviour in the angles implicit in @ and 
the linear dependence with slope 2.211 ± 0.014 (stat.) ± 0.053 (syst.) well 
consistent with the expected slope Ca/Cf HZ]- The data are sensitive to 
the small term oc l/iV^ which corresponds to a negative interference not 
accessible through purely probabilistic jet algorithms. It is remarkable that 
the perturbative calculations describe production of particles with such low 
momenta (few 100 MeV above threshold) at low multiplicity. 

Similar transverse effects are expected for "fp collisions (qq antenna in 
direct, gg in resolved processes) and pp collisions (low transverse radiation 
in direct 7 or W production, high radiation in gluon jet production). 

Underlying event in pp collisions (R.D. Field): Soft particle produc- 
tion in central rapidity perpendicular to the high px trigger jet direction 
increases rapidly with the jet transverse momentum from low px (like mini- 
mum bias) up to about 7 GeV and saturates beyond. The same is observed 
for the transverse momentum sum in back-to-back jet events. If one triggers 
in addition on a particle in transverse direction one observes the "birth" of 
a third jet in the same and possibly even on a fourth jet in the opposite 
direction. This conclusion is derived from the good agreement with the 
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PYTHIA MC which includes multiple parton collisions whereas HERWIG 
without this addition is in a less good agreement. It will be interesting 
in future studies to clarify the role of multi-parton scattering and also to 
investigate the possible reduction of the transverse particle production in 
direct 7 and W production processes as expected for the perturbative gluon 
radiation mechanism emphasized in the previous paragraph for e + e~ colli- 
sions. 



4. Small x structure functions and diffraction 

There is an old expectation for the parton density at small x to "satu- 
rate" [2S], i.e. to become so high that a limiting behaviour related to the 
finite proton (or nuclear) size is reached. Ultimately, one expects a tran- 
sition into a strong coupling regime not accessible to perturbative treat- 
ment. With the new data from HERA and RHIC this debate enters a new 
round. The "Pomeron" which describes diffractive processes with vacuum 
quantum number exchange is treated as composite object with a partonic 
sub-structure. 



4-1. Deep Inelastic Scattering and parton saturation 

In the standard perturbative treatment of DIS (DGLAP 1972-1977) the 
photon interacts with the proton through exchange of a single parton ladder 
which leads to a linear evolution equation of the parton densities in Q 2 . 
With decreasing Bjorken x there is the possibility of multiple interaction 
of photon and proton through exchange of two or more parton ladders as 
described by GLR |28| (1983). This happens for sufficiently large parton 
overlap probability W(x, Q 2 ), given by the ratio of the parton parton cross 

section at scale Q 2 , a ~ as q$ > to the mean distance of two partons in the 

proton, Ab 2 ~ ' f° r proton radius R and parton density F(x,Q 2 ) ~ 

xG{x,Q 2 ). While for W <C 1 the DGLAP approximation is appropriate, 
for W ~ a s the non-linear recombination processes set in and, ultimately, 
at W = 1 saturation is reached, i.e. a full overlap of partons in the proton 
which is beyond perturbation theory. This limit defines the characteristic 
saturation scale Q s (x) from 

xG{x,Q 2 8 ) Q s {x) 2 

ttR 2 ~ a s (Q 2 ) [ 1 



For small a s this corresponds to a state of high gluon density. 
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The theoretical analysis starts from the dipole picture |29| . formulated 
in space time, from which the j*p total cross section can be computed as 

4 P L {x,Q 2 ) = f d 2 r f dza dipolc (r,x)\i;l L (f,z,Q 2 )\ 2 (4) 

where ip' 7 is the wave function of the virtual photon splitting into a qq 
dipole, z the photon longitudinal momentum fraction of the quark and r 
the transverse size of the dipole. 

Different approaches are used for the dipole cross section. The scat- 
tering process can be studied in the proton rest frame and one considers 
higher orders to the qq wave function. A non-linear evolution equation for 
the dipole-proton scattering amplitude has been given by Balitsky and by 
Kovchegov [SO]. A complementary approach treats the gluons at small x 
in an effective field theory in a frame with a low momentum photon and 
an energetic proton where the partonic motions in the proton are largely 
frozen in. The state of these high density gluons is also called "Colour Glass 
Condensate" (CGC) : J^. This new kinematic regime of perturbative QCD 
at high density at the border to a non-perturbative confinement regime is 
under intense investigation and is important also for heavy ion collisions 
where the gluonic state at high density appears initially. Status and appli- 
cations of the theory of saturation and the CGC are reviewed by J. Bartels 
and E.G. Ferreiro. 



4-2. Evidence for saturation at HERA? 

According to the above outline one may reach the saturation region 
in DIS either by decreasing x at fixed Q 2 , i.e. by increasing the par- 
ton density or by decreasing Q 2 at fixed x, i.e. by increasing the parton 
transverse "size". In the first case one observes that the structure func- 
tion i<2 can be well fitted by the NLO QCD in the DGLAP approach 
for Q 2 > 2 GeV 2 where for x < 0.01 one finds an x independent slope 
MQ 2 ) = —{dhxF^/ d\nx)Q2 which depends linearly on Q 2 . For smaller 
Q 2 the applicability of the perturbative calculations may be questioned, so 
there is no direct evidence for saturation in the perturbative regime from 
this point of view. 

On the other hand, a new regime appears in the second case when Q 2 
is decreased at fixed x. In this case it is observed that the A slope saturates 
for Q 2 < 1 GeV 2 (talk by E. Elsen). This region of low Q 2 is included in 
the models of saturation which combine perturbative and non-perturbative 
aspects. Some essential features of this approach are contained already in 
the model by Golec-Biernat and Wiisthoff |32|. Here the DIS cross section 
is obtained in the dipole picture (jlj) with a simple ansatz for the dipole 
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cross section (Jdi VO \ c {r 2 Q 2 s {x)) to depend only on the particular combination 
of r and x and with Q 2 (x) ~ x~ x the saturation scale. The cross section is 
calculated perturbatively for small distances (adipoie ~ r2 ) whereas at large 
distances a simple Gaussian form has been adopted with saturation built 
in (odipoie - > Co)- In this way the full Q 2 range becomes accessible in the 
model. 

An important prediction of the model is the geometrical scaling [SB] 
which states that the cross section aj ^{x, Q 2 ) = /(r) depends only on the 
quantity r = Q 2 /Q 2 (x). This scaling property is well satisfied. 1 Another 
successful prediction of the model is the near constancy of the ratio of the 
diffractive over the full cross section F^ /F^ * under variation of the total 
hadronic mass W . 

Recent studies have removed some shortcomings of the model, especially 
the phenomenological parametrizations for large dipole sizes r. The large Q 2 
behaviour can be recovered by a smooth matching to the DGLAP evolution 
|35| . The solutions of the BK equation determine the behaviour for large r 
and together with the DGLAP behaviour at small r a good description of 
HERA data is obtained : the geometrical scaling in the saturation 

domain, the transition between the hard and soft photon region for the slope 
A and the x-dependence of the saturation scale Q 2 . These successes in the 
description of data involving the scale Q s can be taken as indirect evidence 
for the onset of saturation. 



4-3. Hard diffraction in DIS 

Events with a large rapidity gap adjacent to the proton have been ob- 
served in NC events with high Q 2 at HERA which are considered as inelastic 
diffraction of the photon [K. Borras; C. Kiesling). Similar events with the 
charged current have now been reported as well. Assuming Regge factor- 
ization for small momentum transfer t between the in and outgoing protons 
this process can be described by Pomeron exchange where the virtual pho- 
ton scatters off the Pomeron which is emitted by the incoming proton with 
momentum fraction xjp = M x /s where Mx denotes the hadronc mass of 
the 7*iP and ^fs the cms energy of the j*p system. In analogy to the 
standard parton model for ep DIS one can introduce Parton Distribution 
Functions for the Pomeron [SHI to describe elP DIS. In this "Diffractive 
DIS" it is possible to derive QCD factorization of the photoabsorption cross 



An alternative scaling law has been derived within the framework of a generalized 
vector-dominance model I34| . 
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section 7* p — ► pX at fixed x ip and momentum transfer t 

d 2 a DDIS /-ajp 



dxipdt 



which further simplifies according to Regge factorization f^(xip, t;x,Q 2 ) = 

fiP/p( x iP> t)fq/ip((3 = x/xip, Q 2 )- These PDF's are then studied as function 
of j3 = x q / IP or f3 = Xgjjp where Bjorken x = j3xip. Both factorization 
properties are found to be satisfied by the data but it is necessary to include 
Reggeon exchange in addition to Pomeron exchange. 

The very precise data from the full HERA-I analysis show the rise of 
the reduced cross section a 1 ^((3, Q 2 ) (suitably normalized to correspond to 
F% p ) with Q 2 for (3 < 0.7 and the decrease for higher (3. The change in slope 
occurs at much higher value than in case of j*p scattering. The striking 
pattern of scaling violation is explained in a NLO DGLAP fit (parameters 
Ajjij, initial PDF's at Qq = 3 GeV) by the large gluon fraction in the 
Pomeron f g = 75±15%. QCD factorization is also verified by the observation 
of diffractive di-jet and charm jet production in agreement with NLO QCD 
computations using the Pomeron PDF's so obtained. 

Theoretical models for the diffractive structure functions (J. B artels) 
have been developed within the dipole picture with saturation, such as the 
models by GBW ISS. already mentioned and by BEKW jlU] which take the 
higher order QCD processes (7* — ► qqg) into acount and provide a good 
description of F 2 DD/5 (/3, Q 2 ). 

The structure of simple Feynman diagrams for DDIS has been discussed 
by S. Brodsky. Explicit calculations in case of Feynman gauge for 7*^ — > 
ssq show that the rescattering of the struck s quark involving nearly on- 
shell intermediate states leads to an imaginary amplitude and an effective 
Pomeron exchange in the production of the colour singlet ss state. This 
process survives in the Bjorken limit. The same QCD final state interaction 
also can produce single spin asymmetries in semi-inclusive DIS |41j . 



4-4- Diffraction in pp collisions 

Results from the TEVATRON on multi-gap events, hard diffractive 
processes with jets, W, Z, J ftp, B have been discussed; exclusive double 
Pomeron Xc 77 and di-jet production are of interest as benchmark for 
exclusive Higgs production at the LHC; the analysis of RUN II data is in 
progress (reports by K. Borras, M. Convey, K. Goulianos). 

The factorization of Pomeron processes has been established in DDIS 
but the arguments cannot be taken over to hadronic processes. In fact, 
the observed cross sections for diffractive di-jet rates at the Tevatron, are 
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Fig. 3. NLO cross sections for diffractive dijet photoproduction as functions of jet 
energy fraction x J 7 ots compared to preliminary HI data. The predictions including 
absorption (R=0.34 03]) agree with the data, those without absorption (R=l) do 
not (from Ref. 

suppressed by an order of magnitude as expected using the Pomeron PDF's 
from HERA assuming factorization |42| . 

Such a suppression has been expected from reinteraction of spectator 
partons which yields a reduced gap survival probability, see for example 
|43| . This idea is supported by the observation of CDF that the occurence of 
an additional gap is unsuppressed; by appropriate combination of multi-gap 
cross sections in soft diffraction the survival probability has been determined 
as S = 0.23 ± 0.07 at y/s = 1800 GeV. Taking this effect into account the 
factorization properties and agreement with extrapolations from HERA can 
be reestablished. A general systematics of multi-gap events of various kinds 
with characteristic factorization properties, also in hard processes, has been 
proposed (K. Goulianos) in the framework of a parton model which includes 
empirical rules like "1/M 2 scaling" and "Pomeron flux renormalization" . 

The theoretical approach by KKMR 44 is based on hard QCD scat- 
tering processes, but includes initial state interaction by multiple Pomeron 
exchanges which are derived in a 2-channel eikonal model. This approach 
explains quantitatively the phenomena of factorization breakdown and the 
rates of multiple gap events. An interesting prediction concerns the break- 
down of factorization in di-jet photoproduction at HERA ,441 with an addi- 
tional suppression S = 0.34. This effect has been verified recently in a NLO 
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QCD calculation [l^l in comparison with HI data [SHI (see Fig. EJ). 

5. Heavy ion collsions and Evidence for Quark Gluon Plasma 

Here the transition from a parton to a hadron ensemble is studied in a 
very high particle multiplicity environment which suggests a thermodynamic 
treatment. In lattice QCD one expects with increasing temperature T or 
energy density e a phase transition from confined to deconfined matter, i.e. 
from a hadron gas to a Quark Gluon Plasma. The ratio e/T 4 , a measure 
of the number of degrees of freedom, shows a characteristic rise with T 
near To ~ 170 MeV, eo ~ 0.7 GeV/fm 3 , depending also on the number 
of flavours, over a range of about AT ~ 80 MeV and then, above e ~ 2 
GeV/fm 3 , becomes nearly T independent; in this region there is still a large 
(~ 30%) deviation from the asymptotic Stefan-Boltzmann limit for the ideal 
quark-gluon gas [17]. Another prediction concerns the dependence of the 
critical temperature T c on the baryochemical potential [i which can be tested 
through the hadron composition of the final state. 

Previous research at the SPS has identified various signatures expected 
for the transition to a QGP, such as strangeness excess, J/V'-suppression, 
universal chemical freeze out; the energy density is found at e ~ 2 — 4 
GeV/fm 3 , just above the critical value. Now with RHIC a new regime with 
the much higher initial density of ~ 15 GeV/fm 3 is reached much above the 
critical density. This leads to new signatures: a strong asymmetric flow of 
particles reflecting the initial spacial anisotropy and the strong absorption 
of high px jets in the nucleus ("jet quenching"). It appears particularly 
impressive with the new RHIC data, that more specific QCD tests are now 
becoming feasible. A general outline of the RHIC results and their inter- 
pretation is given by R. Seto. 

5. 1 . Onset of deconfinement in the SPS energy range 

Results on AA collisions from an energy scan over the lower SPS energies 
20-80 AGeV have been presented by P. Seyboth. After the observation of 
signatures for QGP formation at the top SPS energy the aim was to search 
for an energy threshold of such signatures. 

A remarkable effect is seen in the energy dependence of the freeze out 
temperature, as determined from the slope of the transverse mass of kaons: 
at low energies there is a strong rise of temperature followed by saturation 
over the SPS energy range and continued rise at RHIC energies. This is 
the typical behaviour expected for a mixed phase of QGP and hadron gas 
at constant temperature. Such a behaviour at the energy density e > 2 
GeV/fm 3 in the SPS range matches the values expected from lattice calcu- 
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lations. A threshold effect is also seen for strangeness production, especially 
the K + /tt + ratio which is a well established signature for QGP formation. 

Another test of lattice QCD calculations concerns the phase diagram in 
the variables T vs. baryochemical potential [is- The analysis of the particle 
species abundances in statistical models yields values for T and hb at freeze 
out which converge for low \xb (high energies) towards the QCD expectation 

5.2. Space-time evolution of collision process at RHIC 

Initial stage (E.G. Ferreiro, K. Tuchin). The initial conditions can be 
introduced by the gluon density in the nucleus at the saturation scale (see 

Eq. ©) as xG A {x,Q 2 s ) ~ ^J^' A) ("Colour Glass Condensate") where 

the A dependence is inferred from Ga ~ A, nR\ ~ At and Q 2 S ~ As. 
Assuming proportionality of particle multiplicity and initial gluon rapid- 
ity density alN/dy ~ xG(x,Q 2 s ) one finds for an AA collision with N par t 
participating nucleons a very slow increase of central hadron multiplicity 
(1/ Np ar t) dN/dy ~ ^/ a s{Ql) ~ hiAp art with energy and centrality [3H|, a 
very successful prediction; other predictions from parton saturation follow 
for the transverse momentum distributions (review ;50 ). 

Early interactions: jet production and jet quenching (Miller, Vitev). 
These new phenomenona are related to the hard parton-parton scatter- 
ing in nuclear collisions and subsequent absorption of one parton in the 
dense medium. The absorption is mainly due to induced gluon radiative 
energy loss in multiple scattering inside the nucleus and is proportional 
to the plasma density jSJ. The most striking evidence for jet quenching 
comes from the study of azimuthal angle correlations of particles associated 
with a high p? trigger particle [52] : whereas in both pp and dAu scattering 
one observes a jet in direction opposite to the trigger jet this away side jet 
is fully suppressed in central AuAu collisions. This is naturally explained 
by the hard collision taking place near the edge of the nucleus where one 
scattered parton leaves the nucleus undisturbed whereas the second parton 
has to move through the big nucleus. This measurement demonstrates the 
big difference between cold nuclear matter traversed in dA collisions and 
the matter created in the central AuAu collision. Theoretical calculations 
lead to an estimate of energy density of e ~ 15 GeV/fm 3 corresponding to 
about 100 times nuclear density. Additional studies confirm the absorption 
strength as function of the nuclear thickness in non-central collisions as well 
as the reappearence of the lost energy of the primary parton in the soft 
particles in jet direction. 

While these observations provide already a strong argument in favour 
of the presence of a QGP there are further crucial tests. QCD predicts the 
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Fig. 4. Transverse momentum dependence of elliptic flow V2 measured by STAR 
|55j a) before and b) after rescaling by the number of valence quarks n. 

absorpion of a gluon to be stronger than that of the quark by the factor 
Ca/Cf- Furthermore, heavy quarks are less absorbed as the small angle 
radiation is cut off below C = rriQ / E ( "dead cone effect" ) [HSl • This effect 
is being searched for in nuclear D meson production [Z. Xu). 

Hydrodynamic evolution, flow phenomena (Y. Rama, T. Hirano, H. 
Long, J. Velkowska and S. Voloshin). The hydrodynamics description of 
the particle production process includes the initial condition (energy den- 
sity, initial flow), the EoS for the QGP with a phase transition built in 
(parameter: latent heat) and the freeze out mechanism for hadron produc- 
tion. The observations related to hydrodynamics are: 

1. Violation of mr-scaling in AA collisions, which denotes a universal slope 
(3 in dN/dp\ ~ e -P m T f or particles of different mass and works well in pp 
scattering. In hydrodynamic flow the particles acquire similar velocities and 
therefore protons obtain higher momenta than pions. 

2. An important consequence of hydrodynamics is the appearence of an 
asymmetric, especially elliptic flow: in the non-central collision of two nu- 
clei there is an almond shape overlap region which generates an asymmetric 
pressure gradient with maximum in impact direction; this results in a cor- 
responding asymmetry in energy and particle flow ^¥ ~ 1 + 2v2 cos 2(j) + . . . 
The elliptic flow V2 increases for particles of higher px with a delay for heav- 
ier particles. 

3. The Equation of State (EoS) can be investigated, especially the prop- 
erties of the phase transition (latent heat ~ 800 MeV/fm 3 ) which provides 
another QCD test. 

4. An apparent problem for the hydrodynamic description is met with the 
Bose-Einstein correlations between identical particles which depend on the 
evolution of the space time volume containing the matter. Whereas T. Hi- 
rano notes a failure of this description it has been pointed out by T. Csorgo 
in the discussion that the disagreement can be avoided by a proper choice 
of the initial transverse flow to explain the final "Hubble flow" (Buda-Lund 
model pU). 
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Coalescence, quark recombination (R. Hwa, J. Velkovska, S. Voloshin). 
Another striking new phenomenon observed at RHIC is the grouping of 
spectra according to the number of constituent quarks. This is observed in 
the pt dependence of the elliptic flow parameter V2 {jpt < 6 GeV) where 
7T, K and p, A, H fall into separate bands but show a uniform dependence if 
rescaled according to the number ny of valence quarks [SSI 023 

v 2 /nv = fipr/nv)- (6) 

As an example the K, A spectra are shown in Fig. |1J Another observation 
concerns the ratios Rcp of particle spectra for central and peripheral col- 
lisions where mesons (4>, K° , K^) and baryons (fi, E,A+ antiparticles) fall 
into separate bands in the region 2 < px < 6 GeV. This confirms the idea 
of parton coalescence |57| . There are several other "anomalies" in nuclear 
production, for example, the large ratio p/ir + for large pt > 2 GeV which 
can be explained within a recombination mechanism for thermal and shower 
partons (R. Hwa). 

The behaviour (JUJ) suggests that before hadron formation there is a flow 
of constituent quarks which then recombine into the observed hadrons. This 
implies a strong dependence of the q/g composition of the plasma during the 
evolution as illustrated in Fig. [S] : initially the primary collision produces 
mainly gluons at high temperature; during the expansion qq pairs are pro- 
duced, but in approaching the critical temperature the gluons are absorbed 
by the quarks ("constituent quarks") which then by coalescence form the 
final state hadrons. The strong ordering according to valence content is 
against expectations from a hadronic resonance gas which would group par- 
ticles according to mass. These observations are therefore another strong 
argument for the presence of a QGP (although at the end without gluons). 
Looked at in reverse order the evolution depicted in Fig. |S] is quite natu- 
ral: the hadrons under increased pressure dissolve at first into constituent 
quarks (as in the additive quark model) but under increased pressure gluons 
are easily freed from the constituent quarks and yield a genuine QGP. 

Strongly interacting QGP. There is another interesting message in the 
PT-dependence of the elliptic flow parameter v<2- Using a calculation for 
a QGP within transport theory |58| it is found that the gluon-gluon cross 
sections of few mb expected in pQCD would give negligeable flow effects, 
only cross sections of > 40 mb would yield the observed asymmetric flow. 
Therefore, the data suggest a strongly interacting QGP and this can be 
related to the large deviation from the Stefan-Boltzmann limit of the ideal 
gas found in the lattice QCD calculations |59| . 



18 




high T 
gluon 
dominated 



T decreasing 
quarks 
produced 



quark 
dominated 




hadrons 



Fig. 5. Evolution of the quark gluon plasma from a gluon dominated phase of 
high temperature and density towards a quark dominated phase near the critical 
temperature T c , finally the transition from constituent quarks to hadrons. 



Finally, we list a few other topics which have been discussed. 

1. Hadronic Phenomena. Particle correlations are often not accessible to a 
QCD description, especially the Bose-Einstein correlations which turn out 
to be particularly important in the discussion of Heavy Ion Collisions. Also 
there are discussions of critical phenomena such as percolation and cluster- 
ing in a string or hadron model framework. 

2. Hadron Spectroscopy. This is another active field with surprising results 
on unexpected hadronic particles, including the still controversial "pen- 
taquarks", especially 6 + (1540), not observed at TEVATRON, and the new 
charmonium state with decay X(3872) — > J/^7r + 7r~ reported here by CDF, 
which stimulates discussions of nonperturbative aspects of QCD. 

3. Astroparticle Physics. In this field we witness a flourishing activity 
where Multiparticle Production plays an important role, although not as 
the basic goal of the activity but rather as a tool. Primary cosmic parti- 
cles of superhigh energies (10 s TeV) are studied. The interpretation of the 
particle yields, in particular the determination of the primary particle en- 
ergies requires a detailed understanding of the propagation and showering 
of the cosmic rays (elementary particles or nuclei) in various media (air, 
earth, ice), including the saturation phenomena and production and decay 
of heavy quarks (c, b quarks) . Questions of the role of particles in theories 
beyond the standard model are being discussed, as well as the origin of such 
high energy cosmic rays. 



The number of multiparticle production phenomena which can be ex- 
plained within QCD is steadily increasing. 



6. Other presentations 



7. Summary 
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1. Hard processes: Calculations for larger kinematical ranges, for observ- 
ables of higher complexity and with increasing accuracy agree with the data, 
no definitive failures have been reported this time. 

2. Soft particle production: It follows perturbative QCD expectations sur- 
prisingly well which is in support of a parton hadron duality picture of 
hadronization with soft colour confinement. 

3. Small x and diffraction: For a high density regime at small Q 2 there is 
indirect evidence for saturation ("Color Glass Condensate") inferred from 
the success of saturation models. Intrinsic Pomeron structure is a useful 
concept in diffractive scattering, the systematics of factorization and its 
breaking are becoming better understood. 

4- Heavy Ions and QGP: There is an indication of a phase transition with 
a mixed phase over the SPS energy range. The higher initial pressure and 
longer evolution time available at RHIC have provided clear evidence for 
jet quenching and strong elliptic flow with an initial energy density about 
100 times higher than in nuclear matter and an order of magnitude above 
critical density. These phenomena are adequately described in terms of a 
strongly interacting QGP; these observables can be used as new diagnos- 
tic tools which allow detailed tests of (perturbative and non-perturbative) 
QCD predictions: parton type dependence of absorpion, Equation of State 
with latent heat and strong deviation from ideal gas limit, phase diagram 
T c vs. fiB, the initial CGC state. Hadrons are formed by coalescence of 
constituent quarks which dominate the QGP in its final stage. 
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